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The behaviour of turbulent flows within the single-layer quasi-geostrophic (Charney- 
Hasegawa-Mima) model is shown to be strongly dependent on the Rossby deformation 
wavenumber A (or free-surface elasticity). Herein, we derive a bound on the inverse energy 
transfer, specifically on the growth rate d£/dt of the characteristic length scale £ repre- 
senting the energy centroid. It is found that d£/dt < 2 ||<7|| /(£ S X 2 ), where WqW^ is the 
supremum of the potential vorticity and £ s represents the potential enstrophy centroid of 
the reservoir, both invariant. This result implies that in the potential energy dominated 
regime (t > £ s ^> A -1 ), the inverse energy transfer is strongly impeded, in the sense that 
under the usual time scale no significant transfer of energy to larger scales occurs. The 
physical implication is that the elasticity of the free surface impedes turbulent energy 
transfer in wavenumber space, effectively rendering large-scale vortices long-lived and in- 
active. Results from numerical simulations of forccd-dissipative turbulence confirm this 
prediction. 



1. Introduction 

The atmosphere and oceans are hugely influenced by the background planetary rotation 
and stable density stratification. These features virtually suppress motion in the vertical 
direction, rendering geophysical fluid systems, intrinsically vast in their horizontal extent 
compared with their limited depth, approximately two-dimensional (2D) at large scales. 
Many models derived from the theory of quasi-geostrophy exploit this simplification. One 
such model, governed by the Charney-Hasegawa-Mima (CHM) equation, is the subject 
of the present study. 

The CHM equation, which describes the evolution of potential vorticity in a shallow 
fluid layer with a free surface, takes the form (Pedlosky 1987) 

|f + J(V,«) = o. (i.i) 

Here tp(x, y, t), the streamfunction, is assumed to be periodic in both x and y. Note that 
ip is proportional to the free surface height anomaly. The quantity q = (A — \ 2 )ip is the 
potential vorticity, f where A is the Rossby deformation wavenumber (A = f/c, where 
/ is the Coriolis frequency and c is the small-scale gravity waves speed, both assumed 
constant). The differential operators </(•, •) and A are, respectively, the Jacobian and two- 
dimensional Laplacian. Eq. (|1.1|) expresses material conservation of potential vorticity. 

f This form of potential vorticity ignores the so-called /3-term. For some studies concerning 
the effects of this term, see Kukhakin & Orszag (1996), Okunu & Masuda (2003) and Smith 
(2004). 
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It also governs the evolution of quasi-2D fluctuations of the electrostatic potential in a 
plane perpendicular to a strong magnetic field applied uniformly to a plasma, in which 
case ifi is the electrostatic potential and A -1 is the ion Larmor radius (Hasegawa & Mima 
1978). 

The CHM equation has been a subject of active research for decades (see for ex- 
ample, Hasegawa & Mima 1978; Fyfe & Montgomery 1979; Yanase & Yamada 1984; 
Larichev & McWilliams 1991; Ottaviani & Krommes 1992; Kukharkin, Orszag & Yakhot 
1995; Iwayama, Shepherd & Watanabe 2002; Arbic & Flierl 2003; Tran & Bowman 
2003). A majority of these studies have been concerned with the dual transfer of the 
total energy (-ipq)/2 = (|V^| 2 }/2 + X 2 (ip 2 )/2 and potential enstrophy (Aijjq)/2 = 
(|A?/>| 2 }/2 + A 2 (|V?/>| 2 )/2, which are conserved by the CHM nonlinearities. Here (■) de- 
notes a spatial average. The kinetic energy and the usual enstrophy are (\Vip\ 2 )/2 and 
(\Aip\ 2 )/2, respectively (these are not conserved). The quantity A 2 (?/> 2 )/2 is the potential 
energy. In addition to the total energy and potential enstrophy, the CHM dynamics also 
conserve an infinite class of potential vorticity norms, including the L°° norm H^H^- Simi- 
lar to 2D Navier-Stokes (NS) turbulence, the dual conservation law implies that the total 
energy (potential enstrophy) is preferentially transferred to lower (higher) wavenumbers. 
Unlike 2D NS turbulence, each of the invariants in CHM turbulence consists of two com- 
ponents. The presence of the second component is due to the elasticity of the free surface, 
which is known to have a fundamental influence on the dynamics. A remarkable effect is 
that the inverse energy transfer for wavenumbers k <C A (the so-called potential energy 
regime)f is predominantly a transfer of the potential energy (Tran & Bowman 2003), i.e. 
of the quadratic quantity (ip 2 ). 

Numerical studies of CHM turbulence appear to suggest that in the potential en- 
ergy regime, the inverse cascade is strongly impeded. For turbulence decaying from an 
initial energy reservoir, Larichev & McWilliams (1991) observe a "slowness" in the tur- 
bulent evolution. Iwayama et al. (2002) notice the formation of a sharp spectral peak, 
a manifestation of an inverse energy transfer with limited extent in wavenumber space. 



Arbic & Flierl (2003) find that the inverse cascade slows down and that vortex merg- 



ing takes a very long time. For forced turbulence, Kukharkin et al. (1995) observe the 
formation of vortical "quasicrystals" , states of long-lived inactive vortices. While these 
observations are consistent with the finding of Tran & Bowman (2003) that the inverse 
energy transfer involves virtually no kinetic energy when k <C A, they are not fully ex- 
plained by this fact. A more satisfactory explanation would be that the inverse transfer 
of potential energy is also impeded in this regime. 

In this Note we derive an upper bound for the growth rate of a characteristic length 
scale associated with the inverse energy transfer. We infer from this upper bound that 
in the potential energy regime, the inverse energy transfer is strongly impeded by the 
elasticity of the free surface. An interesting physical interpretation of this effect is that 
this elasticity prohibits the merger of large-scale vortices to form larger energy-carrying 
scales, effectively rendering such vortices long-lived and inactive. Although this effect is 
deduced in the framework of unforced and inviscid dynamics, it also manifests itself in a 
forced-dissipative version of (|1.1|) . as is illustrated by some numerical results. 



| For scales smaller than A -1 = c//, the motion is inherently three-dimensional (cf. Dritschel 
& de la Torre Juarez 1996; Dritschel, de la Torre Juarez & Ambaum 1999) and therefore poorly 
described by l|l.l[l . 
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2. Inverse energy transfer 

The nonlinear transfer of energy to lower wavenumbers of an energy reservoir (when 
it occurs) necessarily entails an increase of a suitably defined length scale representing 
the energy centroid of the reservoir. Here we consider the length scale I defined by 

(H) (|v^| 2 ) + A 2 (v^ 2 } ' 1 ■ 1 

In H2.1|) the operator (— A) e is defined by (— A) e ip{k) = k 2e i}j{k), where k — \k\ is the 
wavenumber and ip(k) the Fourier transform of tp(x). The growth rate d£/dt provides 
a quantitative measure of the inverse transfer. In particular, a smaller rate corresponds 
to an inverse transfer with a more limited extent in wavenumber space than one with 
a greater rate. In this respect, d£/dt may be considered as a measure of locality of the 
inverse transfer, in the sense that for a given energy distribution a smaller rate implies a 
more local transfer than a greater one. 

In the case where most of the kinetic energy resides in wavenumbers k <C A (the 
potential energy regime), (|(— A) 1 / 4 -;/;! 2 } is negligible compared with A 2 (|(— A) -1/,4 ?/>| 2 ), 
and similarly (|V-!/>| 2 ) is negligible compared with X 2 (ip 2 ). More precisely, in this regime 
we have 

{K-A)vyi 2 ) (|v^| 2 ) 

A 2 (|(-A)- 1 /V|2) - x 2 (^ 2 ) ' 

where the first inequality is due to Holder's inequality (see also Tran & Shepherd 2002 and 
Tran 2004), implying £ w (|(— A)~ 1 / 4 ip\ 2 )/(ijj 2 }. For an energy reservoir in the potential 
energy regime, the length scale £ specifies where, in wavenumber space, most of (ip 2 ) (or 
equivalently most of the potential energy) is distributed. Its growth rate is thus the rate 
at which the spectral profile of the potential energy proceeds toward lower wavenumbers. 

The evolution equation for the characteristic length scale £, defined in the preceding 
paragraphs, is obtained by multiplying l|l.lf> by (—A)~ 1 / 2 ip/(ipq) and taking the spatial 
average of the resulting equation: 

ld£ _ {(-A)~ 1 / 2 tpJ(ip, q)) (gJQM-A)- 1 / 2 ^)) 
2dt (|W>| 2 ) + A 2 (^ 2 ) (|V^| 2 ) + A 2 (^ 2 ) 

<|g||V^||V(-A)-V 2 ^|) (1^2)1/2(^1/2 

(|V^| 2 ) +A 2 (V 2 ) " (|V^| 2 ) + X 2 (yj 2 ) 1191100 ' K '' } 

where the various manipulations of the nonlinear term are straightforward. In particular, 
the last step can be readily seen by expressing the streamfunction in terms of a Fourier 
series. 

The ratio in front of ||g|| on the right-hand side of l|2.2|l can be closely estimated in 
terms of conserved quantities. For that purpose, let us denote by £ s the conserved length 
scale representing the enstrophy centroid of the reservoir: 



H/>9> 1/2 _ «|V^| 2 ) + A 2 ^ 2 }) 



1/2 



(A^q) 1 / 2 {{\AtP\ 2 ) + A 2 (|VV>| 2 )) 1/2 ' 
Solving (|2.3|l for (ijj 2 ) and substituting the result into the said ratio we obtain 

Qvyf) 1 / 2 ^ 2 ) 1 / 2 _ (iV^^^dA^ + ^A 2 -!)^! 2 )] 1 ^ 
(|V^| 2 ) + A 2 (^ 2 } A [(| Wf) + ^ 2 (|AVf) + {£ 2 X 2 - 1)(|V^| 2 )] 



(2.3) 
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(\viA 2 ) 1/2 ((\^\ 2 ) + * 2 m\ 2 )) 1/2 

4A((|A^|2) + A 2 (|Wf}) 

im 2 ) l/2 



4A((|A^|2) + A 2 (|Wf}) 



< Tv>- (2-4) 



V2 " £\2 



This is a very sharp estimate, especially for a reservoir in the potential energy regime. 
We note that in general the inequality I > i s holds. This is a consequence of the Holder 
inequality and can be shown as follows. Let us denote by U(k) the power spectrum of 
{— (— A) -1 / 2 ^). The length sales t and i s can then be expressed in terms of U(k) as 

._ JU(k)dk 



and 

is = 

Now by Holder's inequality we have 



JkU(k)dk 

J kU{k) dk 
Jk 3 U(k) dk 



1/2 



/ r \ 2/3 / f N 1/3 

kU(k)dk<( U{k)dk\ I k 3 U(k)dk 

from which the inequality £ > £ s indeed follows immediately. 
Putting these results together we have 

f t < 2 JMk < 2 |!k, (2 . 5) 

using £ > £ s . Hence the exponential growth rate of £ is bounded from above by 2 WqW^ / [i 2 A 2 ). 
Given a finite WqW^, the upper bound 2 WqW^ /(^ 2 A 2 ) can be made arbitrarily small pro- 
vided ^ 2 A 2 is sufficiently large. This condition is satisfied for an initial energy reservoir 
in a sufficiently remote wavenumber region of the potential energy regime, where both 
£ S X ^> 1 and £X ^> 1. Under these circumstances, no significant growth of £ occurs when 
the initial energy reservoir freely spreads out, and hence no significant inverse energy 
transfer is possible. In passing it is worth mentioning that the upper bound for the expo- 
nential growth rate of £ in 12.5fl consists of the constant A and the two invariants WqW^ 
and £ Sl and therefore can be completely determined by initial data. 

The result in the preceding paragraph is an analytic version of the qualitative statement 
in the literature that the turbulence evolves on a rescaled (slower) time (cf. Larichev & 
McWilliams 1991; Kukharkin, Orszag & Yakhot 1995). To clarify this point let us rewrite 
(j21>|) as 

df 

Ps*\ t < 2 Moo*- (2.6) 

The left-hand side of (|2.6|l is the (instantaneous) growth rate of £ on the slow time r = 
t/(£ 2 s \ 2 ). For an initial reservoir with a finite I^H^, the inverse transfer of potential energy 
can be considerable only on this slow time scale. In the limit £ S X — * oo, the invariants 
become {tp 2 )/2 and (|Vi/'| 2 )/2 and the length scales £ and l s become (K-A)" 1 / 4 ?/'! 2 }/^ 2 } 
and {ip 2 } 1 / 2 / dVipl 2 } 1 / 2 , respectively. The rescaled time t becomes infinitely long. 

An estimate which may be significantly sharper than H2.5(l can be derived by replacing 
the potential vorticity q in l|2.2|l by the relative vorticity Atjj. By doing that, instead of 
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12.511 . we obtain 



In remote regions of the potential energy regime, where £X > £ S A 3> 1, we expect A ip to 
dominate Aip. Therefore in these regions it is plausible that WqW^ 3> HAt/?^, allowing 
1)2.7(1 to become significantly sharper than 1)2.5(1 . The price for this improvement is that 
II^V'lloo i s n0 ^ conserved. Note that in such regions t and £ s are the length scales of 
potential and kinetic energy, respectively. 

A straightforward physical interpretation of the above result is that the elasticity of 
the free surface suppresses energy transfer for scales larger than the Rossby deformation 
radius. This effect has a profound influence on the vortex dynamics. Let us consider a hy- 
pothetical flow consisting of vortices (having smooth potential vorticity — no small length 
scales) of comparable scales, which are much larger than the deformation radius A . For 
these scales the inverse energy transfer is strongly impeded, allowing for virtually no for- 
mation of larger energy-carrying scales. This keeps the vortices from merging, effectively 
establishing a state of long-lived inactive vortices. This may explain the crystallisation 
of vortices observed by Kukharkin et al. (1995). The impeded inverse energy transfer 
implies a similar effect on the direct transfer of potential enstrophy, a consequence of 
the dual conservation law. One can see from the "symmetry" of the dual conservation 
law that the available potential enstrophy for direct transfer depends on the inversely- 
transferred energy. Similarly, the available energy for inverse transfer depends on the 
directly-transferred potential enstrophy. Hence impeded inverse energy transfer means 
impeded direct potential enstrophy transfer (and vice versa). This implies that smooth 
potential vorticity with scales ^> A -1 remains smooth. 



3. Numerical results 

We now consider results obtained from numerical simulations of a forced-dissipative 
version of ((1.1|) . These convincingly illustrate the strong reduction of the inverse transfer 
in the potential energy regime, and may help explain an observation by Iwayama et 
al. (2002) in their numerical simulations of freely decaying CHM turbulence wherein 
an initial energy reservoir evolves into a single sharp peak. A plausible explanation for 
this observation is that the sharp peak is due to both the viscous dissipation (which 
suppresses the formation of small scales) and the impeded inverse transfer deduced in 
the preceding section: their collective effect "squeezes" energy into a narrow band of 
wavenumbers, resulting in the observed peak. In our forced-dissipative case we observe 
that in the potential energy regime the inverse energy flux is arrested and deposits most 
of the potential energy at the lower-wavenumber end of the quasi-steady energy range. 
We also observe a relaxation of the kinetic energy spectrum of the inertial range (from 
fc -5 / 3 toward fc -1 ) as the inverse cascade proceeds from k > A to k < A. 

We simulate Eq. ((1.1(1 . with the addition of molecular viscosity and a body force, in a 
doubly periodic square 2n x 2tt with a dealiased 1024 2 pseudospectral method. No large- 
scale dissipation mechanisms are employed and no steady dynamics are sought since we 
are primarily interested in the near arrest of the inverse energy cascade in a quasi-steady 
state. (For an analysis of steady CHM turbulence in the presence of large-scale damping, 
see Smith et al. 2002.) Three different values of A are used: A = 20, 40 and 60. The forcing 
is spectrally localised in the wavenumber interval K = [49.5,50.5], in the sense that its 
Fourier components /(fc) are nonzero only for those modes k having magnitudes k lying 
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Figure 1. The kinetic energy spectrum E(k) vs. k averaged over periods [139, 140], [155, 156] 
and [160, 161] for A = 20, 40 and 60, respectively. 



in the interval K: 

/( fc ) = ± . (3.1) 

N E \^p)\ 2 
i P i=fc 

Here e = 1 is the constant energy injection rate and N is the number of distinct wavenum- 
bers in K. The viscosity coefficient is v — 5 x 10~ 4 . The simulations were initialised with 
the kinetic energy spectrum E(k) = 10 _5 7rfc/(50 2 + k 2 ). 

Figure ^ shows the quasi-steady kinetic energy spectrum E{k) vs. k averaged over 
periods [139, 140], [155, 156] and [160, 161] for A = 20, 40 and 60, respectively. In all three 
cases, the accumulation of energy at low wavenumbers is clearly visible. This is attributed 
to the impeding effect discussed above, and not finite-size effects since the energy levels 
at the lowest wavenumbers are still small compared with the energy peaks. This is true 
for both kinetic and potential components. For A = 20 we have a limited kinetic energy 
regime in the wavenumber interval [20,50]. Before the energy peak reaches A = 20, 
the dynamics of this kinetic energy regime are expected to resemble those of 2D NS 
turbulence: the energy cascades to low wavenumbers via a fc -5 / 3 kinetic energy spectrum. 
We indeed observe a discernible fc -5 / 3 spectrum (not shown) in the wavenumber range 
[20, 50] before the energy peak reaches A = 20. The dynamics after the energy peak 
passes A undergo profound changes. Most notably, the growth of potential energy occurs 
mainly around the energy peak, consistent with the reduction of the inverse energy 
transfer in the conservative case. Another remarkable change is in the spectral slope of 
the energy spectrum from k = A = 20 up to the forcing wavenumber k = 50. The slope 
of the kinetic energy spectrum in this region appears to relax from —5/3 to — 1. The 
quasi-steadyt spectrum becomes fc _1 with an excessive accumulation of energy at the 
lower-wavenumber end. This unexpected feature is common to the other two cases, for 



f Note that after the energy peak passes A, growth of energy is due primarily to that of the 
potential energy and the kinetic energy becomes quasi-steady. 
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Figure 2. The energy centroid I versus time t (left) and scaled time i(20/A) 2 (right) for 

A = 20, 40 and 60, as labelled. 

which hardly any inverse transfer region exists between the forcing wavenumber and A 
and for which no clear fc -5 / 3 spectrum is formed initially. 

The evolution of the energy centroid I for these three cases is shown in Figure [21 on 
the left versus unsealed time t and on the right versus scaled time i(20/A) 2 . Both the 
reduction in d£/dt with increasing A and the collapse of the curves when plotted versus 
scaled time are consistent with the growth rate bound i|2.5[l .t Only the A = 20 case stands 
out slightly, due to an early kinetic energy transfer which is virtually absent in the other 
two cases. 

4. Summary and discussion 

The main result of this work is an upper bound for the growth rate of the characteristic 
length scale associated with the inverse energy transfer in CHM turbulence. This upper 
bound is expressible in terms of the supremum of the potential vorticity, the characteristic 
length scale representing the potential enstrophy centroid of the reservoir and the Rossby 
deformation radius. These are constants of motion, allowing the bound to be completely 
described in terms of initial data. It is found that in remote regions of the potential energy 
regime, the turbulent transfer from an initial energy reservoir becomes strongly impeded, 
in the sense that hardly any potential energy is transferred to lower wavenumbers. This 
effect has been illustrated by numerical simulations of a forced-dissipative version of the 
CHM equation. The physical implication of the present finding is that the elasticity of 
the free surface hinders the nonlinear transfer of potential energy to larger scales. 

An interesting implication of the present result is that vortices having smooth potential 
vorticity and scales much larger than the Rossby deformation radius are kept from merg- 
ing to form larger energy-carrying scales and remain smooth. Such vortices are inactive 
and long-lived compared to their counterparts in 2D NS turbulence. 

The present theory is fully consistent with other results for both forced (Kukharkin et 
al. 1995) and freely decaying (Larichev & McWilliams 1991; Iwayama et al. 2002; Arbic 
& Flierl 2003) CHM turbulence, as discussed above. These, we argue, are explained by 
the growth rate bound (|2.rJ|l and its physical interpretations given at the end of §2. 

Analyses of atmospheric and oceanic data yield results that could plausibly be ex- 

f We cannot directly compare the numerical growth rates with the theoretical bound because 
the imposed forcing violates conservation of WqW^ and l a used in 
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plained by the present work. Boer & Shepherd (1983) find that the spectral flux in at- 
mospheric wind data is arrested with essentially no inertial range. |Scott fc Wang (2005)| 



find the inverse energy transfer in the oceans to be qualitatively similar to the result of 



Boer & Shepherd (1983) These are in remarkable agreement with the idealised numerical 
results noted above and with the present findings. 

Iwayama et at (2002) observe "well-behaved" tails of the potential vorticity probability 
density function (PDF) — much steeper than observed in NS turbulence, see their figure 
5 — and raise the question why it is that CHM turbulence in the potential energy regime 
exhibits no intermittency. A speculative answer provided by the authors is based on an 
asymptotic version of the CHM equation wherein the relative vorticity is dropped from 
the time-derivative term: 

A 2 ^ + J(A^) = 0. 

This equation, as noted by Iwayama et al. (2002), describes advection oiip by Aip (i.e. by 
the smaller scales of motion). These authors argue that this may lead to more "random 
walk" behaviour, effectively suppressing intermittency. The present study provides an 
alternative answer: the reduction of turbulent transfer both renders large-scale vortices 
inactive and discourages erratic small-scale motion. This also can account for the observed 
steep PDF. In some sense CHM turbulence in the potential energy regime is far less 
turbulent than its 2D NS counterpart. 

The numerical simulations in § 3 were carried out, using a computer code developed by 
John Bowman, when CVT was a Pacific Institute for the Mathematical Sciences postdoc- 
toral fellow at the University of Alberta. Constructive comments from two anonymous 
referees were very much appreciated. 
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